UNITARITY AND INTERFERING RESONANCES IN 

7T7T SCATTERING AND IN PION PRODUCTION ttN -> ttttN. 

M. Svec* 
Physics Department, Dawson College, Montreal, Quebec, Canada H3Z 1A4 

and 
Physics Department, McGill University, Montreal, Quebec, Canada H3Z 1A4 

o 

(N 

£ 

3 



OO 



x 



Abstract 
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(— 1 , of Breit-Wigner amplitudes with specific complex coefficients which depend 

' on the resonance parameters of all contributing resonances. We use analyt- 
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icity of 7T7T partial wave amplitudes to show that they must have the form of 



a coherent sum of Breit-Wigner amplitudes with complex coefficients and a 
complex coherent background. The assumption of additivity of Breit-Wigner 
phases is a new constraint on partial wave amplitudes independent of par- 
tial wave unitarity. It restricts the partial waves to analytical functions with 
very specific form of residues of Breit-Wigner poles. Since there is no phys- 
ical reason for such a restriction, we argue that the general form provided 
by the analyticity is more appropriate in fits to data to determine resonance 
parameters. The partial wave unitarity can be imposed using the modern 
methods of constrained optimization. We discuss the production amplitudes 
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in ttN — ► 7T7riV reactions and use analyticity in the dipion mass variable to 
justify the common practice of writing the production amplitudes in pro- 
duction processes as a coherent sum of Breit-Wigner amplitudes with free 
complex coefficients and a complex coherent background in fits to mass spec- 
tra with interfering resonances. The unitarity constraints on tttt partial wave 
amplitudes with resonances determined from fits to mass spectra of produc- 
tion amplitudes measured in irN — ► ttttN reactions can be satisfied with an 
appropriate choice of complex residues of contributing Breit-Wigner poles. 
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I. INTRODUCTION 

In 1930's, Breit and Wigner introduced MM a parametrization of resonances observed in 
the energy dependence of integrated and differential cross-sections of nuclear reactions. The 
original Breit- Wigner formula was only a one-resonance approximation and its justification 
was initially only phenomenological. A theoretical justification for Breit- Wigner formula 
later emerged from quantum collision theory ||. The evident existence of multiple and 
overlapping resonances in nuclear reactions led to two distinct generalizations of the Breit- 
Wigner formula for an isolated resonance to multiresonance description of the scattering 
process. 

One generalization was undertaken by Feshbach HH, Humblet || and McVoy @ who 
used the analyticity properties of the S'-matrix to show that the transition matrix can be 
written as a coherent sum of Breit- Wigner terms with complex coefficients and a coherent 
background. Since the transition matrix must satisfy unitarity, the parameters and coeffi- 
cients of this multiresonance parametrization are not independent [|5]|§. In principle it is 
possible to use the methods of nonlinear programming PJlOf and constrained optimization 



with computer programs such as MINOS developed at Stanford University (TT| to impose 
the conditions of unitarity in fitting the experimental data. 

Another approach to multiresonance description of scattering process was proposed by 
Hu in 1948 ||12| . He observed that the Breit- Wigner contribution of an isolated resonance 



to the S'-matrix is unitary and proposed to describe the multiresonance contributions in the 
S-matrix by the product of isolated Breit- Wigner contributions for each resonance. Since 
each term is unitary, the product also satisfies unitarity. The partial wave phase shift is 
then a sum of Breit- Wigner phases of contributing resonances and a background phase. 
As a result, the expressions for partial wave amplitudes involve products of Breit- Wigner 
amplitudes. This method has been recently used by Bugg et al |L3| and by Ishida et al |L4j 
in their analyses of ttti phase shift data. 

Up to now the connection between these two descriptions of multiresonance contributions 



(interfering resonances) has not been clarified. In this work we show that the Hu description 
is a special case of a more general description based on analyticity. We show that the Hu 
method also leads to a coherent sum of Breit-Wigner amplitudes with complex coefficients 
and a complex coherent background for any partial wave as expected from the analyticity 
of the S-matrix. However, the complex coefficients have a very specific form in terms of 
resonance parameters of all contributing resonances. The assumption of additivity of Breit- 
Wigner phases is a new constraint that restricts the partial waves to analytical functions 
with these specific residues of Breit-Wigner poles. Furthermore we show that the additivity 
of Breit-Wigner phases is an assumption entirely independent of the unitarity property of 
partial wave amplitudes which is a condition imposed on their inelasticity. 

Since there is no physical reason why the physical partial waves must have the form of 
a coherent sum of Breit-Wigner amplitudes with the specific complex coefficients required 
by the additivity of Breit-Wigner phases, we conclude that the general form imposed by 
the analyticity is more appropriate for fits to data to determine resonance parameters. This 
conclusion is particularly relevent for analysis of interfering resonances in the mass spectra in 
production processes such as 7rN — > ttttN or pp — > irirpp. Using analyticity in the invariant 
mass variables we justify the common practice of parametrizing the production amplitudes 
in terms of a coherent sum Breit-Wigner amplitudes with free complex coefficients and a 



complex coherent background 15-23 



The paper is organized as follows. In Section II we briefly review the unitarity and 
the problem of interfering resonances in potential scattering since it motivates the analysis 
in hadronic reactions. In Section III we review the two-body partial wave unitarity in tttt 
scattering and its relation to the general form of isospin partial waves. In Section IV, we 
introduce the assumption of additivity of Breit-Wigner phases in the tttc scattering and 
show that it leads to partial waves in a form of a coherent sum of Breit-Wigner amplitudes 
with specific complex coefficients and a coherent background. In Section V we generalize 
dispersion relations for partial wave amplitudes in tttt scattering to Breit-Wigner poles and 
show that the form obtained from the additivity of Breit-Wigner phases is a special case. 



In Section VI we focus the discussion of the two methods to a finite energy interval and 
argue that the addition of Breit-Wigner phases imposes an unjustified constraint on fits to 
data. In Section VII we formulate unitarity for production amplitudes in ir~p — > 7T~7T + n 
reaction and contrast it with partial wave unitarity in tttt scattering. In Section VIII we 
show that the method of addition of Breit-Wigner phases can be generalized to production 
amplitudes. We also use analyticity in the invariant mass to obtain a more general form 
for production amplitudes in terms of a coherent sum of Breit-Wigner amplitudes with free 
complex coefficients (pole residues) and a complex coherent background. We argue that this 
general form is more appropriate in fits to measured mass spectra. Although the discussion 
is confined to pion production amplitudes in ttN — > ttttN, the conclusions have general 
validity. We also comment on determination of tttt partial wave amplitudes from resonance 
parameters determined in measurements of production amplitudes in ttN — > ttttN reactions. 
The paper closes with a summary in Section IX. 

II. UNITARITY AND INTERFERING RESONANCES IN POTENTIAL 

SCATTERING. 

A. Unitarity 

We will consider the scattering of a spinless particle of mass m by a real, central potential 



V(r) [p4J1 - In the asymptotic form of the stationary scattering wave function, the outgoing 
wave is characterized by the scattering amplitude f(k, 9) where k is the wave number of the 
particle related to its energy by 

E = ^ = — (2.1) 

2m 2m y ' 

and 9 is the scattering angle. In the units h — 1 the wave number k has the meaning of 
momentum p. The scattering amplitude can be written in the form 

oo 

f(k 1 0) = y £{2e+l)T t (k)P t (co8e) (2.2) 



The partial wave amplitudes Tg are given by 

Tl = ^k [Si{k) ~ 1] (2 ' 3) 

where Sg(k) is called ^-matrix. For elastic scattering 

S t = e 2l5i{k) (2.4) 

where the phase-shifts Se describe the interaction and are related to the potential V(r). For 
elastic scattering \Se\ = 1 which is the condition of elastic unitarity. 

When a particle collides with a target, non-elastic processes are possible and particles 
are removed from the incident (elastic) channel. Since the interaction can alter only the 
outgoing part of the wave function, we require that the amplitude of the outgoing wave 
be reduced if non-elastic processes occur. The reduction of scattering amplitudes leads to 
conditions of inelastic unitarity 

\S e \ < 1 (2.5) 

This suggests that we write 

S e = m (k)e 2iS ^ (2.6) 

where r\e is called inelasticity and has values 

< w < 1 (2-7) 

The partial wave then has a general form 



;From (2.8) it follows that 



I m Tz = k\T t \ 2 + ±-{l-4) (2.9) 



4A- 
This equation expresses the unitarity condition on the partial waves T t 



B. Interfering Resonances. 

In the following we will work with partial wave amplitudes 

1 



U = m = ^[ m e 2i5i -l} (2.10) 



and the energy E instead of k. A detailed study of the potential scattering [23] shows that 
the phase shift may be decomposed as 6t = £e + pe where ^ is the background phase which 
does not depend on the shape and depth of the interaction potential V(r) while the part p? 
does depend on the details of the potential. Near resonant energy E r 

tm ^2(!^ (2 - U) 

where T(E) is the width of the resonance. We introduce a Breit-Wigner resonance phase 5 r t 

51 = tan ~ 1{ 2(E r ( ~V = arg[Er ~ E+ 4 r(jB)] (2 - i2) 

such that in the energy interval AE centered about E r we have pi ~ b\ and 

SiaZi + Si (2.13) 

Prom (2.12) it follows that 

„2* S r __E-E r -\iT _ 1 + 2 . ai (214) 



E- E r + \iY 



where 



a '=E^TW (2 - 15) 

is the Breit-Wigner amplitude of the resonance r. For an isolated resonance we then obtain 

it = J- fa ^- 1)+ ^ ( _g£ > (2 , 6) 



If N resonances contribute over an interval AE then, following Hu ]nj and references |p~3| , pT4 
we can write 



e mi = tt *v a' 1 = TT (1 + 2ia n) ) (2.17) 

i 1 ! £ - E^ + liTW n=l 

The prescription (2.17) clearly satisfies unitarity but seems to lead to a complicated expres- 
sion for partial waves tn in terms of Breit-Wigner amplitudes af . 

On the other hand, analyticity of S- matrix was used by Feshbach |3]|| , Humblet |j) and 
McVoy [0] to derive a general form for tt || 

t e = B t {E) + V H-^ 1 ( 2 - 18 ) 

tiE-E^ + \iT^\E) 

where Bf is a background term and A™ (E) are complex coefficients. The sum in (2.18) can 
be written as a coherent sum of Breit-Wigner amplitudes 



IV 



i r W(E) 



t t (E)=B t (E) + y Rf\E) ^ — ^ (2.19) 

K ' K ' ti J E-Ei n) + ^TM(E) V ' 

In Section IV we show that the prescription (2.17) leads to the analytical form (2.19) with 

specific expressions for the coefficients Rp' and the background Bg. 



III. ISOSPIN AMPLITUDES AND UNITARITY IN vrvr -+ vrvr SCATTERING. 

Hadron resonances have definite values of spin and isospin. It is therefore necessary to 
express the amplitudes for charged pion processes tttt — > 7r7r in terms of isospin amplitudes 
T ! (E, 9) with definite isospin / = 0, 1, 2 and work with partial wave amplitudes T[(E) p5| . 
At first we will work with the center-of-mass energy E = \/i to pursue the analogy with the 
potential scattering. 

The partial wave amplitudes T/ satisfy partial wave unitarity equations |25| J/ 



ImTl = q\Tl\ 2 + Al (3.1) 

where A{ are the contributions from inelastic channels, such as 7tty — > nnnn, KK, NN, and 



q = J\(s — 4/i 2 ) is cm. momentum where li is the pion mass.. Let us write A\ in the form 

A, 7 = ^(l-(^) 2 ) (3.2) 

8 



Then the unitarity equation (3.1) has the same form as (2.9) and the partial waves T/ can 
be written as 

Tl = ^Wi*** 1 ~ 1} (3.3) 

where the 5{{E) are phase shifts and the inelasticity 



4{E) = Vl - 4qAi(E!) (3.4) 

is given by the inelastic unitarity contributions A|. In analogy with potential scattering we 
expect that < rjj < 1. As we shall see later, the descriptions of interfering resonances in 
7T7T scattering do not depend on the condition that r]\ <\. 

The positivity of inelasticity rjj in (3.4) imposes a constraint 

A,' < 1 (3.5) 

We can now show that the unitarity equation (3.1) admits no solution for A| > l/4g. If the 
inelastic unitarity contributions satisfy this condition we can write 

A^ = l(l + (^) 2 ) (3.6) 

where (r/j) 2 > 0. Setting T/ = V/ /2q we get from the unitarity equation (3.1) 

(ReV/) 2 + (ImV/ - l) 2 = -(ri) 2 (3.7) 

which is not possible. Thus the conditions (3.5) represent genuine constraints on the inelastic 
unitarity contributions Aj and the parameterization (3.3) of partial wave amplitudes with 
(3.4) is the most general solution of the unitarity equation (3.1) for all E. 

Since the values of inelastic terms A| in the partial wave unitarity equations (3.1) are 
not known, we constrain the partial waves T/ by inequalities imposed by the unitarity. From 
the positivity of r}\ and the condition (3.5) we obtain 

ImTl < q\T\\ 2 + 1 (3.8) 



If we add the requirement that r\\ < 1, then A\ > from (3.2), and we also have the usual 
unitarity constraint 

ImTl > g|T/| 2 (3.9) 

The inequality (3.9) implies positivity 

ImT* > (3.10) 

at all energies. 

Finally we note the following observation. Let f(z) be any complex function. Then 
1 + 2if(z) is a complex function that can be written as 

1 + 2if(z) = r](z)e 2lS ^ (3.11) 

where f] > and S is real. Thus any complex function can be written in the form 

f(z) = \\j)e 2i8 - 1) (3.12) 



and satisfies the equation 



Imf=\tf + \(l-r?) (3.13) 



We see that the unitarity equations (3.1) are a special case of (3.13) with rj given by (3.4). 

IV. INTERFERING RESONANCES IN vrvr SCATTERING USING THE 
ADDITION OF BREIT-WIGNER PHASES. 

The general form of phase shift parametrization of partial wave amplitudes T/ is 

T/ = i-[^e 2 ^-l] (4.1) 

with inelasticity r\\ determined by unitarity via (3.4). In the following we will omit the 
indices t and / for simplicity. In analogy with potential scattering, we decompose the phase 
shifts 5 into two parts 
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5 = $ + S r (4.2) 

where £ is the nonresonant background phase and <5 r is the phase due to physical particle 
resonances occurring in the partial wave T/. The phase of a single isolated resonance is 
given by the Breit-Wigner formula 



2i8 r _ E - E r - \iT(E) 



(4.3) 



"" E - E r + \iV{E) 

Let us consider that N resonances contribute to the partial wave amplitude T/. Following 
ref. [12-14] we now assume, that the resonant phase shift 5 r is given by the sum of the 
Breit-Wigner phases of the contributing resonances 



N 



£ ^ (4-4) 



n=l 

We assume that N is finite. Then 



N N 



e 



E-E„- kT 



= n e 2 < = n e-eJ* (4 - 5) 

n=l n=\ ^ - C/ ™ ' 2° n 

We can write for each Breit-Wigner phase 

e 2< = 1 + 2 ian ( 4 _g) 

where a n is the Breit-Wigner amplitude 



2 l r 



e — E n + 2«r r 

Then we can write 



(4.7) 



e 2iS- = x + 2^^ ( 4 _g^ 

where T res is given in terms of products of Breit-Wigner amplitudes a n The partial wave 
amplitude T[ then has a general form 

T=^-( V e 2 ^-l) + -e^ V T res (4.9) 

Let us consider the case N = 2. Then the resonant part of the amplitude T/ is 

11 



T res = Oi + a 2 + 2iaia 2 (4-10) 

where the interference term 

(-i ri )(-|r 2 ) 

0102 " (E-E 1 + \iT 1 )(E-E 2 + ^ir 2 ) (4U) 

With a notation 

z fc = £ fc - -*r fc (4.12) 



we write 



4 +lA-& (4.13) 



(E- Zl )(E-z 2 ) l E- Zl E-z 2 'C 
The requirement that this equality holds leads to relation 

(A + B)E - (Az 2 + Bzx) = C (4.14) 

Next we require that A = —B to eliminate the E dependent term, and get C = A(z,\ — z 2 ). 
Then (4.13) has the form of a sum 

1 11 1 

(E-z l )(E-z 2 ) ~ z 1 -z 2 }E-z 1 ~ E-zJ ( ' 15) 

and we can write the resonant part (4.10) of the partial wave amplitude as the sum of two 
Breit-Wigner amplitudes 

Tres = C { 2) ~^\ + Cf ] =t? (4.16) 

1 E - E 1 + hT x 2 E-E 2 + tT 2 K } 



where the complex coefficients 



«1 ~ Z 2 

are exactly such that the unitarity condition 

12 



C{ 2) = \-2i^^- (4.17) 

z\ - z 2 



C? = 1 + <li -$± 



J2i& r 



|l + 2iT res | = l (4.18) 



is satisfied for all E. The energy dependence of the widths T n (E) introduces energy depen- 
dence in C n (E),n =1,2. 

Consider now the case of three interfering resonances N = 3. Then 

T res = oi + a 2 + a 3 + 2i(aia 2 + a\a 3 + a 2 a 3 ) + (2i) 2 aia 2 a 3 (4.19) 

We can write the last term as a sum 

1 A B C l 

(E - Zl )(E - z 2 )(E - z 3 ) ~ [ E- Zl + E-z 2 + E-zJl5 ( ^ 

Requiring that the terms proportional to E 2 and E in the numerator on r.h.s. of (4.20) 
vanish, we obtain a sum 

Kl + j^ + j<^_ (421) 



where 



E — Z\ E — z 2 E — z 3 



A z 3 — z 2 , 

K i = p = -^Y^ ( 4 -22) 



T s B z 3 - zi 
K2 = D = ~1T 

C z\ - z 2 



K 3 = ^ 



D X 



with 



X = ziz 2 (zi - z 2 ) + z 3 z x {z 3 - zi) + z 3 z 2 (z 3 - z 2 ) (4.23) 

The resonant part of the partial wave amplitude is again a coherent sum of the Breit-Wigner 
terms with complex coefficients 

^ = £Cf ~} V \. V (4.24) 

„=1 E- h n + 2 ?i n 

where 
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Cj 3 ) = 1 - 2i^- - 2i^- + {2i) 2 {\y 2 ){\y z )K 1 (4.25) 

Z\ — Z2 Z\ — Z% Z Z 

cf = i + 2* -£l- _ 2 * -JEl_ + ( 2i ) 2 (ir 1 )(ir 3 )K 2 

Zl — Z2 Z2 — Z% Z Z 

C< 3) = 1 + 2*^- + 2i -ti- + (2z) 2 (ir 1 )(ir 2 )^ 3 

This procedure is general and valid for any finite N. Assuming that the resonant phase 5 r 
can be separated from the phase shift S and is given by the sum of Breit-Wigner phases, we 
will always get the resonant part T res of the partial wave amplitudes T/ in (4.9) as a sum of 
Breit-Wigner amplitudes 

r ~w=£™^?bj (4 - 26) 

In (4.26) the complex coefficients C^ have an explicit form in terms of resonance parameters 
E n ,T n ,n = 1,...,N such that T res satisfies the unitary condition (4.18). The form of 
coefficients C^ depends on the number of resonances contributing to the partial wave T/. 
As the result of (4.26) we can conclude that the multiresonance parametrization of partial 
wave amplitudes based on additivity of Breit-Wigner phases has a general form of a coher- 
ent sum of Breit-Wigner amplitudes a n with complex coefficients and a complex coherent 
background 

1 N -±F 

T = t[B(E) + £ Rl N \E) F 2 " r ] (4.27) 



where 



fl = Ifo e *€-l] (4.28) 



RW = e**nCW = (1 + 2»B)Cf) 

Comparing (4.27) with expression (2.19), we see that the description of multiresonance 
contributions using the addition of Breit-Wigner phases leads to the same form of partial 

14 



wave amplitudes as the analyticity of the S'-matrix. However, the complex background B 
and the complex coefficients R^ in (4.27) have the explicit form (4.28) imposed by the 
additivity of Breit-Wigner phases. 

Note that in the derivation of (4.26) for T res , and in the resultant form (4.27) with (4.28), 
we have not needed or used the assumption that inelasticity i] < 1. The Hu method is based 
on the unitarity of 1 + 2iT res and is not related to the unitarity of the whole partial wave 
amplitude T/. 

Finally we give a relativistic form for the multiresonance description of partial wave 
amplitudes. The relativistic form of Breit-Wigner amplitudes (4.7) is given by 

-m n YJs) 



s 2 -ml + im n T n (s) 



(4.29) 



where we have used m n instead of E n to emphasize that E n is the mass of the resonance. 
To obtain the corresponding coefficients C^\ we make replacements in (4.17) or (4.25) 

-r n -»• m n T n (4.30) 

1 

Z n J-j n 7) n ^"n ''"n ifiTl n L n 

The partial wave amplitudes then have the relativistic form 

1 N 

T(s) = -[B(s) + Y,R { n N) (s)a n ( S )] (4.31) 

where B and -R„ -* are still given by (4.28) with replacements (4.30) to satisfy the unitarity 

V. GENERALIZED DISPERSION RELATIONS FOR PARTIAL WAVE 
AMPLITUDES AND INTERFERING RESONANCES IN tttt SCATTERING. 

In this section we shall relate the multiresonance parametrization (4.31) of partial wave 
amplitudes T/ with a multiresonance parametrization obtained from analyticity. To this 
end we shall use generalized dispersion relations for the amplitudes 

15 



t^s) = qTl{s) (5.1) 

where s is the Mandelstam energy variable. 

Our starting point is the well-known |27| dispersion respresentation of a complex function 
f(z) with simple poles at z n ,n = 1,2, ... ,N in the complex plane z, a branch cut along a 
positive real axis from a to oo and with asymptotic property |z|/(,2) — > as \z\ — > oo. We 
shall also assume that the function f(z) is a real function f(z*) = f*(z). Using Cauchy's 
integral theorem and the process of contour deformation, it can be shown |27j that 



/(*) = f -^- + i IM^M (5.2) 

±< Z - Z n IT J x' - Z 

IL—± a 

A remarkable feature of the proof of (5.2) is that it takes place for a fixed value of z [[27|]. As 
the result, the dispersion relation (5.2) is also valid for moving poles for which z n = z n (z). 
In such a case the residues R n in (5.2) also depend on z, i.e. R n = R n (z). Furthermore, 
the dispersion relation (5.2) is easily generalized to include a left-hand cut and for functions 
that are not real. In the latter case Imf(x') in (5.2) is replaced by a discontinuity function 
along the cut(s). 

In 7T7T scattering, the partial wave amplitudes t\(s) have a right-hand cut for s > 4/z 2 
(where // is the pion mass), and a left-hand cut for s < due to Mandelstam analyticity 
p8| . Let us assume that the amplitude t\ has a finite number N{ of complex poles 



m 



im n T(s), n = 1, . . . ,Nl (5.3) 



corresponding to the resonances in t\. Note that the imaginary part of the poles depends 
on the energy variable s. In principle, the mass m n could also depend on the energy s. This 
possibility has been recently considered by Pennington p5[ . Omitting the indices / and £, 
the generalized dispersion relations for the partial wave amplitude t\ read 

t(s) = I(s) + E -^- (5.4) 



where I(s) are the dispersion integrals over the left-hand and right-hand cuts [2^| and 
R n (s) are the pole residues. It is convenient to rewrite (5.4) in a form using Breit-Wigner 
amplitudes a n (s) 

16 



N 

t(s) = I(s) + J £R n (s)a n {s) (5.5) 

n=l 



where we have redefined the pole residues with 

( \ ~ m n*- n{ s ) /-. a \ 

a n{s) = y— ^y^ 5.6 

s -m 2 n + im n T n (s) 

The representation (5.5) is valid for all s > 4{i 2 . The representation (5.5) of partial waves 

t\ coincides with the parametrization (4.31) provided that 

I(s) = B(s) = ±[r,(8)e*M - 1] (5.7) 



R n (8)=R^(8)= V (8)^^C^( l 



[i + 2iB)a 



(N) 
n 



We see that the multiresonance parameterization based on additivity of Breit-Wigner phases 

(4.4) imposes a special form on the dispersion integrals and pole residues given by (5.7). 

In general, a partial wave t\ can be written in two forms 

1 N 

t = Tplve 215 -l] = I+J2 R n a n (5.8) 

Z% n=l 

Apart from the partial wave unitarity equations (3.1) and (3.4) and the analyticity as- 
sumptions, there are no constraints on the partial waves. The assumption of additivity of 
Breit-Wigner phases (4.4) is a new constraint that restricts the partial waves to analyti- 
cal functions that satisfy the conditions (5.7). We find no physical justification for such 
a restriction and no advantage in using it in phenomenological fits to data to determine 
resonance parameters. 

VI. INTERFERING RESONANCES IN A FINITE ENERGY INTERVAL. 

In the previous two sections we have assumed that N is the total number of resonances 
contributing to a partial wave. The parametrizations (4.31) and (5.5) were valid for all ener- 
gies s > 4/j 2 . In practice iV is not known and fits to data are done in a finite energy interval. 
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Such is the case e.g. of analyses |T3| . |T4}| . In this Section we develop parametrizations of par- 
tial waves appropriate for analyses in a finite energy interval where only a few resonances 
contribute. The parametrizations will be based both on additivity of Breit-Wigner phases 
and analyticity and we shall compare their use in practical fits to data. The results will be 
used in the Section VIII. 

Let us consider energy interval 4/i 2 < s < sm where M resonances contribute. In the 
framework of the assumption of additivity of Breit-Wigner phases we will assume that the 
resonant phases of resonances outside of this energy interval are absorbed in the background 
phase. The total phase shift then is 

5 = i (M )+5 r M (61) 

where 

M 

5li = E 5 : (6-2) 

n=\ 

N 

£<")=£+ J2 % (6.3) 

n=M+l 



The partial wave then takes the form 



i = I(^ , -l)+^ (M) «) (6.4) 



where the resonant part 








T (M) 
res 


M 

= E cm* 

71=1 


is unitary 








e 2iS M 


= i + 2iT%; 



Alternatively we can rewrite (4.31) in the form 



M 



(6.5) 



(6.6) 



t = £?( A/Af ) + ]T R { n N) a n (6.7) 



n=l 



where 

B&* = 1(^-1)+ £ ijgOo™ (6.8) 

/? m=A/+l 

is the background term. Note that the sum in (6.7) is not unitary. We cannot compare the 
coefficients of Breit-Wigner amplitudes a n , n — 1, . . . , M in (6.4) with those in (6.7) since 
£( M ) j n (6.4) contains the terms a m , m = M + 1, . . . , N but the sum in (6.7) does not. 
If we look at the general form (5.5) from analyticity, then for s < sm we can write 

M 

t = £?( M ) + J2 R n a n (6.9) 

ra=l 

where the background term 

N 

BM = I+ £ Rmam (6.10) 

m=M+\ 

In (6.9) the residues R n are not constrained by the conditions (5.7). 

In fitting data using the parametrization (6.4) we explicitly make use of the assumption 
of additivity of Breit-Wigner phases. This is also the case when we use (6.7) if iV is known 
and the coefficients C^' in R^' can be calculated. In general N is not known, and the 
background l?( MiV ) and residues R^ are free parameters. Then there is no difference in 
using (6.7) or the general form (6.9) from analyticity alone, since in (6.9) the background 
£j( M ) anc } residues R n are not constrained except for unitarity. In all cases we use constrained 
optimization of the \ 2 function. In the case of (6.4) we require that inelasticity function 
7] < 1. In the case of (6.7) or (6.9) we require that Imt > |t| 2 and use programs such as 



MINOS JTT| for constrained optimization. 

It is not obvious that the use of parametrization (6.4) from additivity of Breit-Wigner 
phases and the parametrization (6.9) from analyticity alone, will lead to the same resonance 
parameters in both cases. The use of parametrization (6.4) confers no phenomenological 
or computational advantage over the parametrization (6.9). The assumption of additivity 
of Breit-Wigner phases restricts the background and the complex coefficients multiplying 
the Breit-Wigner amplitudes a n , n = 1, . . . , M in the parametrization (6.4) to specific 
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forms. Since there is no physical justification for such a restriction and the parametrization 
(6.9) is free from such constraints, we suggest that the use of parametrization (6.9) is more 
appropriate in determining resonance parameters in mr scattering. 

VII. UNITARITY IN PION PRODUCTION n~p -» vr vr+n 

It is a common misconception to identify the partial wave production amplitudes in 
reaction 7i~p — > 7r~7r + n with partial waves Tj in nn scattering and demand that the partial 
wave production amplitudes also satisfy the partial wave unitarity (3.1). In this Section 
we clarify the distinction between the two kinds of amplitudes and the associated unitarity 
relations. 

The production process n~p — > n~7i + n is described by production amplitudes p5| , |30|j3T 



T XnfiXp (s,t,m 2 ,9,0) (7.1) 

where X p and A n are proton and neutron helicities, s is the c.m.s. energy squared, t is the 
momentum transfer between the incident pion and the dipion system (7r~7r + ), m? is the 
dipion mass squared, and Q = (6, <fi) is the solid angle of the final n~ pion in the dipion 
rest frame. The dipion state does not have a definite spin. The production amplitudes 
(7.1) can be expressed in terms of partial wave production amplitudes M xx ox (s,t,m 2 ) 



corresponding to definite dimeson spin J using the angular expansion [Ea,pQ. 31 



CO +J 

Tx nfi x p = E E (2j + iy/ 2 M^ ni0Xp (s,t,m 2 )di (ey X4, (7.2) 

j=0X=-J 

where J is the spin and A the helicity of the (it~7t + ) dimeson system. 

It is evident from (7.2) that the partial wave production amplitudes M XXn 0A (s, t, m 2 ) 
cannot be identified with the tttc partial wave amplitudes Tj(m 2 ). The amplitudes 
M xx 0A (s, t, m 2 ) can be thought of as two-body helicity amplitudes for a process 7r~ +p — > 
M(J,m) + n where the "particle" M(J,m) has spin J and mass m. 



The production amplitudes T Xnj0Xp satisfy the unitarity condition |26 
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- ^(T A „,oA p - T* Xp< J = £ / T 0Ap , n T; nin d$ n (7.3) 

where (i$ n is the n-body Lorentz invariant phase space of the intermediate state n. Since 
the initial state in 7r~p —* 7r~7r + n is a two-body state and the final state is a three-body 
state, the amplitude T\ nt o\ p enters the unitarity integral only linearly. This occurs only 
when the intermediate state is ir~p or 7r~7r + n. However the three-body intermediate state 
involves a 3 — > 3 amplitude and a three-body phase space integral. Separating the two-body 
intermediate states ir~p and n°n, we can write (7.3) in the form 

- i(Tx nMp - T* Xp>Xn ) = W T 0V0A , T; n>0A , rf$ 2 + W ro Ap , A^; niOA , rf$ 2 + A A „, 0Ap (7.4) 

V J \l J 



\ 



where T 0Apj0A ' and T OXO y n are helicity amplitudes of reactions n~p — > ix~p and n~p — > 7i°n, 
respectively. The amplitude T x ox , corresponds to process n n — > 7i^7i + n. The inelastic 
unitarity contribution A Ani0Ap (s, t,m 2 , 9, (p) can be expanded in the form analogous to (7.2) 

oo +J 

Aa„,oa p = £ £ (2J+l^ 2 Ai XnfiXp (s,t,m 2 )di (e)e^ (7.5) 

j=0\=-J 

Using expansions (7.2) and (7.5) in (7.4) we get unitarity relations for partial wave produc- 
tion amplitudes 

- mL,ox p - M % iXX J = E/ T oA p ,oA p M£ OA ,d$ 2 + J2JT 0XpfiXfn M^ 0K d^ 2 + A{ nfiXp 



K 



(7.6) 



Using time-reversal relations for two-body helicity amplitudes |14] 



ML v ,xx n = (-l) A ^ A ^ A M A J A „, 0Ap (7.7) 

we see that the left hand side of the partial wave unitarity relation (7.6) does not simplify 
to 2ImM xx ox as is the case for the partial waves T/ in tttt scattering. The right hand side 
of (7.6) involves M xx ox only linearly and not quadratically as is the case in inr scattering. 
Futhermore, the right hand side of unitarity relation (7.6) includes (linearly) partial wave 
production amplitudes for the process 7i°n — > 7t~7f + n. We conclude that the unitarity 
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relations (7.6) for partial wave production amplitudes M XXn 0A are complex relations that 
do not have the simple form 

ImTl = q\Tt\ a + Al (7.8) 

of the partial wave unitarity relations in tctc scattering. 

For brevity let us define M/ = M XXn 0A where A stands for the helicities. The amplitude 
M^ is a complex function and so is the function 1 + 2iqMf L . In analogy with (3.12) and 
(3.3) we can write 

Mi = 7^(vie 2iSi ~ 1) (7-9) 

where r]i(s,t,m 2 ) is the "inelasticity" and 5i(s,t,m 2 ) is the "phase shift". The amplitude 
Mj[ satisfies relation similar to (3.13) 

ImMi = q\Mi\ 2 + 1(1 - (r^) 2 ) (7.10) 

Unlike in tttt scattering, the form of the equation (7.10) does not coincide with the form 
of partial wave unitarity (7.6) and the "inelasticity" 7] J k cannot be related to the inelastic 
unitarity contributions A^, in contrast to (3.4). 

VIII. INTERFERING RESONANCES IN PRODUCTION PROCESSES. 

The amplitudes describing the production processes such as ttN — > nnN, pp — > impp 
or pp — > 37r are far more complex than the isospin amplitudes in the tttt scattering. As an 
example, consider pion production in 7t~p —* 7r~7r + n. The angular distribution of the dipion 
7r~7r + state is described by partial wave production amplitudes M( Xn0X (s,t,m 2 ) defined 
in the previous Section with eq. (7.2). The measurements of n~p — > 7T~7T + n on polarized 
target actually determine the moduli of nucleon transversity amplitudes pO|,|3~T| which are 



linear combinations of helicity amplitudes M XXn ox . For masses m <1000 MeV, the pion 
production is described by two >S-wave (J = 0) nucleon transversity amplitudes S and S, 
and by six P-wave (J = 1) nucleon transversity amplitudes L,L,N,N,U,U |l30| , |3~l| . The 
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amplitudes A = S,L,N,U correspond to nucleoli transversity "up" while the amplitudes 
A = S,L,N,U correspond to nucleon transversity "down". The amplitudes L, L correspond 
to dimension helicity A = while N, N and U, U are natural and unnatural parity amplitudes 
corresponding to combinations of A = ±1. 

The measurements of pion production on polarized targets enable to advance hadron 
spectroscopy from the level of spin-averaged cross-sections to the level of spin-dependent 
production amplitudes. These measurements determine the mass spectra \A\ 2 and |v4| 2 of 
spin-dependent production amplitudes. Measurements of rc^p — > 7r~7r + n on transversely 
polarized targets were done at CERN at 17.2 GeV/c [32"H35] and at ITEP at 1.78 GeV/c 



. Measurements of 7T + n — > tt + tt p |[l , 37| -3jj| and K + n — > K + n p |40|j4l|| at 5.98 and 



11.85 GeV/c on transversely polarized deuteron target were also done at CERN. More 
recently it has been shown that mass spectra of production amplitudes can be obtained 
in measurements of n~p — > 7r°7r°n, n~p — > rjrjn [fEJ and 7r~p —> T]7r~n, ir~p — > r]7c°n |43 



on transversely polarized targets allowing for amplitude spectroscopy of these interesting 
processes. 

The analysis of mass spectra measured in production processes requires a parametrization 
of the production amplitudes in terms of the Breit-Wigner amplitudes to identify contribut- 
ing resonances and to determine their parameters. Here we discuss two approaches, one 
based on the additivity of Breit-Wigner phases and the other on analyticity of production 
amplitudes A(s, t,m 2 ) in the mass variable at fixed s and t. 

First we note that the unitarity equation (7.6) for the partial wave production amplitudes 
in ttN — > ttttN is a complex relation and that the helicity amplitudes M/ or the transversity 
amplitudes A and A do not satisfy the two-body partial wave unitarity equation (3.1) with 
(3.4). Nevertheless, the experimentally measured production amplitudes A(s,t,m 2 ) are 
complex functions and as such can be written in the form 

A(s,t,m 2 ) = ^[r ]A e 2 ^-l] (8.1) 

where the "inelasticity" tja = VA(s,t,m 2 ) and "phase shift" 5a = 8A(s,t,m 2 ) depend also 
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on the helicities or transversities of the amplitude A. Obviously, 

I mA =\A\ 2 + ^l- V 2 A ) (8.2) 

However, there is no requirement now that ^ < 1 since tja has no relation to unitarity as 
in (3.4). 

We can pursue the analogy with the im scattering, and impose an assumption that the 
"phase shift" 

8 A (s, t, m 2 ) = U(s, t, m 2 ) + 8 r (m 2 ) (8.3) 

where 5 r is the sum of Breit- Wigner phases of the N resonances contributing to the amplitude 
A and £ A is the "background" phase. If we restrict ourselves to a finite mass interval 
4/x 2 < m? < m 2 M with M resonances, we can write 

8 A (s,t,m 2 ) =& M \s,t,m 2 ) + 6 r M (m 2 ) (8.4) 

A(s, t, m 2 ) = ^(VAe 2 ^ - 1) + ^'^'(m 2 ) (8.5) 

in analogy with (6.4) for mi scattering. 

A more general approach is to use analyticity of A(s, t, m 2 ) in m 2 with s and t fixed. We 
can assume that kinematical singularities have been removed from the production amplitudes 
A(s, t, m 2 ) [45 J . Assuming that there are N Breit- Wigner poles in the amplitude A(s, t, m?) 



in the mass variable m 2 , we can use the generalized dispersion relations for the variable m 2 
with s and t fixed to get 

N 

A(s, t, m 2 ) = I(s, t, m 2 ) + J2 Rn(s, t, m 2 )a n (m 2 ) (8.6) 

n=l 

where I is the contribution of dispersion integrals, R n are complex pole residues, and a n 
are the Breit- Wigner amplitudes (5.6). In a finite mass interval 4/i 2 < m 2 < m 2 M with M 
resonances we can write 

M 

A(s, t, m 2 ) = B (M \s, t, m 2 ) + J2 R n(s, t, m 2 )a n (m 2 ) (8.7) 



n=l 
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where the background 

TV 

B^ M \s,t,m 2 ) = I + ]T R m a m (8.8) 

m=M+l 

We note that for M = N we get back the constraints (5.7) with replacements rj — > tja and 
£ — > £4. Again, the assumption of additivity of Breit-Wigner phases restricts the production 
amplitudes to analytical functions that satisfy the constraints (5.7). 

The measured mass spectra \A\ 2 can now be fitted either with the parametrization (8.5) 
or with the more general parametrization (8.7). There are no unitarity constraints to be 
imposed on the production amplitudes \A\ 2 during the fits since the right hand side of the 
unitarity relation (7.6) is not known and the partial wave unitarity (3.1) or (7.8) for tctt 
amplitudes T/(s) does not apply to the production amplitudes A(s,t,m 2 ). The unitarity 
constraint (3.1) or (7.8) can be imposed only in the analysis of data on the tttt — > tttt reaction 
and below we discuss its effect on tttt amplitudes. 

Since there is no physical justification for the assumption of additivity of Breit-Wigner 
phases in (8.3) and since the form (8.5) confers no phenomenological or computational 
advantage over the more general analytical form (8.7) , we conclude that the use of the form 
(8.7) is more appropriate in fits to mass spectra in production processes to determine the 
resonance parameters of interfering resonances. 

The parametrization of production amplitudes in terms of a coherent sum with complex 
coefficients and a complex coherent background as in (8.7) has been an accepted practice for 
a long time. Such parametrizations first appeared in connection with the possible double- 



pole character of the A 2 meson [15] and the splitting of the Q resonance in K + tt tt + mass 



spectrum [JTJJ. Recently such parametrization has been used in the study of <r(750) — /o(980) 



interference in .S-wave production amplitudes in tt p — > tt tt + u measured on polarized target 



at CERN |T7|)[18| and in the study of a — /o(980) interference in the central collision pp — > 
7r°7r°pp [plj. More recently, an analysis of S'-wave production amplitudes from threshold to 



2 GeV in pp — > 7i°7i°pp was made using three |SIJ and four [[HJ] interfering Breit-Wigner 
amplitudes and a coherent background. The GAMS collaboration used four interfering 
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Breit-Wigner amplitudes and a coherent background in their fit of S-w&ve mass spectrum 
from threshold to 3 GeV in n~p — > 7r°7r°n measured at 100 GeV [E2]. Also recently, the 



Fermilab E791 Collaboration used the form (8.7) to fit Dalitz plot of D + — * tt~tt + tt + decays 
in their search for a scalar resonance a |£J . 

Finally we comment on the determination of tttt partial wave amplitudes from measure- 
ments of ttN — > ttttN. The resonance parameters from the fits to mass spectra such as those 
measured in ttN — > tt + tt~N on polarized targets 0,|l^,^] or in tt~p — > 7i°7i°n |pO|-p^] must 



be the same in tttt partial waves. However, the tttt partial wave amplitudes are expected to 
satisfy the partial wave unitarity constraints (3.1) and (3.4), or rather the inequalities (3.8) 
and (3.9) which for the amplitudes t\ defined in (5.1) read 

\t I e \ 2 <Imt I e <\t I i \ 2 + ~ (8.9) 

The unitarity conditions (8.9) can always be satisfied by an appropriate choice of background 
and complex residues R n (s) in the general parametrization (6.9) based on analyticity Al- 
though the tttt partial waves and production amplitudes in ttN — ► ttttN with the same spin 
and isospin share the same Breit-Wigner poles, they are different analytical functions and 
thus the residues of the poles and the backgrounds are different. In particular, the residues 
in production amplitudes A depend on particle helicities and kinematic variables s and t. 
Accepting the resonance parameters obtained from the fits to the mass spectra \A\ 2 mea- 
sured in ttN — > ttttN to describe the resonances in tttt — > tttt scattering, the effect of the 
unitarity conditions (8.9) is to constrain the residues R n (s) and the background term in the 
general parametrization (6.9) of the tttt amplitudes. 

It is also possible to use the resonance parameters determined from measurements of 
ttN — > 7T7riV to calculate the resonant part T^p and to define the tttt partial waves using 
the parametrization (6.4) with free background and inelasticity functions £} ■ \s) and r](s). 
The unitarity can be satisfied by imposing the condition rj < 1. 

The unitarity constraints (8.9) may not uniquely determine the background and pole 
residues in the parametrization (6.9) from analyticity and the use of parametrization (6.4) 
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from additivity of Breit-Wigner phases is questionable. We conclude that the resonance 
parameters determined from mesurements of ttN — > 7nrN alone may not determine the tttt 
partial wave amplitudes without additional assumptions or direct measurements of tttt — > tttt 
reactions. 

IX. SUMMARY. 

We have shown in the case of tttt scattering that the assumption of additivity of Breit- 
Wigner phases in a partial wave amplitude leads to a sum of Breit-Wigner amplitudes with 
complex coefficients and a coherent background (4.31). The coefficients have a specific form 
(4.28) in terms of resonance parameters of all contributing resonances. The form (4.31) is a 
special case of the general form (5.5) based on analyticity and it is not related to the unitarity 
property of partial waves ((3.1) and (3.4)). The claims |T3| , |T4|| that the additivity of Breit- 



Wigner phases provides the correct description of interfering resonances in tttt scattering 
are not justified since there is no physical reason why the Breit-Wigner poles must have 
the specific residues imposed by this assumption. We found that the Breit-Wigner phases 
of interfering resonances are not necessarily additive. We suggest that the general form 
(5.5) from analyticity is more appropriate in fits to data. The unitarity conditions (8.9) 
m\ 2 < Imt\ < \t\ | 2 + I can be effectively imposed using the modern methods of constrained 
optimization [pHiT|. 

Mass spectra in production processes are described by production amplitudes. We used 
the case of ttN — > tt + tt~N reaction to illustrate the complexity of production amplitudes. 
Specifically, the production amplitudes do not satisfy the two-body partial wave unitarity 
equation (3.1), depend of particle helicities and on several kinematic variables in addition 
to the invariant mass. We have used analyticity of production amplitudes in the invariant 
mass variable to justify the common practice |T5|-p3"f of writing the production amplitudes 
as a coherent sum of Breit-Wigner amplitudes with free complex coefficients and a complex 
coherent background in fits to measured mass spectra to determine the resonance param- 
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eters of interfering resonances. The two-body unitarity constraints on the tttt partial wave 
amplitudes with the same resonances can be satisfied by an appropriate choice of complex 
residues of the contributing Breit-Wigner poles. This reflects the fact that the 7m partial 
wave amplitudes and production amplitudes while sharing the same resonances are different 
analytical functions. 
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